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Abstract. Polarimetric observations of snow surfaces, ob- 
tained in the 41 0-2264 nm range with the Research Scan- 
ning Polarimeter onboard the NASA ER-2 high-altitude air- 
craft, are analyzed and presented. These novel measurements 
are of interest to the remote sensing community because 
the overwhelming brightness of snow plagues aerosol and 
cloud retrievals based on airborne and spacebome total re- 
flection measurements. The spectral signatures of the polar- 
ized reflectance of snow are therefore worthwhile investigat- 
ing in order to provide guidance for the adaptation of algo- 
rithms currently employed for the retrieval of aerosol prop- 
erties over soil and vegetated surfaces. At the same time, 
the increased information content of polarimetric measure- 
ments allows for a meaningful characterization of the snow 
medium. In our case, the grains are modeled as hexagonal 
prisms of variable aspect ratios and microscale roughness, 
yielding retrievals of the grains’ scattering asymmetry pa- 
rameter, shape and size. The results agree with our previous 
findings based on a more limited data set, with the majority 
of retrievals leading to moderately rough crystals of extreme 
aspect ratios, for each scene corresponding to a single value 
of the asymmetry parameter. 


1 Introduction 

The development of accurate techniques for the retrieval 
of climatologically relevant parameters in snow-covered re- 
gions is of obvious importance for the success of climate 
model predictions (Comiso, 2006; Zwally and Giovinetto, 
2011; Zatko et ah, 2013). For instance, the retrieval of grain 
size and surface temperature leads to information that can 
be used to map the snow melting state, valuable to under- 
stand the mass balance of the ice sheets (Stamnes et ah, 


2007; Lyapustin et ah, 2009). Passive optical sensors are 
employed in this regard to evince the optical and micro- 
physical properties of snow grains from measurements of 
snow reflectance, and advanced algorithms have been ap- 
plied to data collected by the MODerate-resolution Imag- 
ing Spectroradiometer (MODIS) and sensors with similar 
capabilities (Hori et ah, 2007; Aoki et ah, 2007; Jin et ah, 
2008; Stamnes et ah, 2011). However, the limited informa- 
tion content of such measurements translates into unique- 
ness problems during the retrieval process, with the result 
that operational snow products are still limited to snow cov- 
ered area and albedo (with the exception of the recent ad- 
dition to MODIS of the dust radiative forcing product as 
described in Painter et ah, 2012). Polarimetric remote sens- 
ing (Diner et ah, 2007; Cairns et ah, 2009; Dubovik et ah, 
2011; Hasekamp et ah, 2011; Chowdhary et ah, 2012; Ot- 
taviani et ah, 2012a; Knobelspiesse et ah, 2012) adds to 
MODIS -type observations because it can provide informa- 
tion on the ice grain shape (van Diedenhoven et ah, 2012; 
Ottaviani et ah, 2012b) and an accurate characterization of 
the optical and microphysical properties of aerosols present 
in the scene, including possible impurities embedded in the 
snowpack. The benefits of passive polarimetric observations 
originate from both their extreme sensitivity to particle mi- 
crophysics and the specifics of surface polarized reflectance. 
For conventional land types (soil, vegetation), the latter can 
be effectively parameterized as a scaled Fresnel reflection 
function with negligible spectral dependence. As a conse- 
quence, the separation of the signal into its surface and atmo- 
spheric contributions (i.e., the atmospheric correction proce- 
dure) is achieved with improved accuracy. 

There certainly is debate on the applicability of the ra- 
diative transfer theory to such a densely packed medium as 
snow. However, the formalism is recognized to model the 
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total reflectance with sufficient accuracy for reasons likely 
associated with the large number of scattering events tak- 
ing place within the snowpack. Based on this assumption, 
the method is expected to perform even better for the po- 
larization component because the polarimetric signatures of 
a medium are known to originate from its top layer (van 
Diedenhoven et al., 2013): deeper into the medium (i.e., past 
the first units of optical depth) multiple scattering random- 
izes polarization. 

The NASA GISS Research Scanning Polarimeter (RSP; 
Cairns et al., 1999), born as an aerosol research instrument, 
provides a target’s total and polarized reflectance with high 
accuracy in nine bands between 410 and 2264 nm and at 152 
viewing angles per scene. Specifically, the dimensionless re- 
flectances are obtained from the measurements of the first 
three components of the Stokes vector (/, Q and U) as fol- 
lows: 
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where the (annual average) extraterrestrial solar irradiance 
Fq is corrected for the Sun-Earth distance ro and the solar 
zenith angle 0 S . 

A recent deployment led to the first published study of the 
polarimetric signatures of snow reflectance from airborne ob- 
servations (Ottaviani et al., 2012b, hereafter referred to as 
012). In that case, the measurements were acquired over 
a topologically complex terrain in the Yosemite National 
Park region and processed with an automated, iterative atmo- 
spheric correction scheme purposely created to analyze the 
spectral dependence of the polarized reflectance. A prelimi- 
nary retrieval of crystal microphysical and optical properties 
was also tested. The purpose of this new work is to apply 
the same atmospheric correction procedure and to extend the 
surface properties retrieval scheme to a more recent data set 
of improved quality, in terms of the flatness of the surveyed 
terrain. 

In the next section we describe the data collection, and 
in Sect. 3 the processing scheme. In Sect. 4 we discuss the 
results obtained from the analysis. The conclusions are found 
in Sect. 5. 


2 Data collection 

The data analyzed here were collected onboard the high- 
altitude ER-2 aircraft during the POlarimeter Definition Ex- 
periment (PODEX) mission, a NASA effort to ultimately 


test and compare the performance of different polarime- 
ters, as part of the Decadal Survey activities and in prepa- 
ration of the Studies of Emissions and Atmospheric Com- 
position, Clouds and Climate Coupling by Regional Surveys 
(SEAC 4 RS) campaign. Compared to more traditional aircraft 
such as the King Air B200 used in 012, the ER-2 offers the 
advantages of a perspective more similar to that of a satel- 
lite. For example, interference caused by cirrus clouds is 
readily assessed using the RSP 1880nm channel, which is 
usually very dark because of the strong absorption by tro- 
pospheric water vapor that masks the surface, and is there- 
fore very sensitive to high-altitude clouds. There is in fact 
interest in exploiting this wavelength as a cirrus screener in 
a broader context, since cold places tend to have dry atmo- 
spheres and the reflectance of snow at 1880nm is much less 
than at, say, 1380nm (the effective strength of absorption by 
ice at 1880 nm is twice that exhibited at 1380 nm). Also, RSP 
operations on the ER-2 benefit from the stable attitude intrin- 
sic to stratospheric flights, which enhances data quality since 
the RSP scans along the flight track (approximately 50° on 
either side of nadir) to assemble the multiangle view of a tar- 
get’s reflectance from subsequent scans. 

About a dozen snow-covered targets were selected for 
the 6.8 h flight over Colorado and the Sierra Nevada range 
in California on 15 March 2012. Well-formed snowpacks 
with little to no interference from vegetation and flatness 
of terrain were the main priorities. Three scenes were cho- 
sen for this analysis. The Grand Mesa plateau (39°08'N, 
108° 12' W, ~ 3000ma.s.l.) in western Colorado exhibits re- 
markable flatness and was covered by a consistent snowpack 
with a depth of greater than 1 m as reported by nearly simul- 
taneous in situ surveys (Mesa Lakes snow telemetry (SNO- 
TEL) site, managed by the Natural Resources and Conserva- 
tion Service of the US Department of Agriculture; M. Skiles, 
personal communication, 2013). Lake Granby (40° 10' N, 
105°53' W, ~ 2500ma.s.l.), the third-largest body of water 
in Colorado and second-largest reservoir, was frozen at the 
time of the overpass and covered with a thin layer of aged, 
likely metamorphosed snow (the SNOTEL site at Stillwater 
Creek recorded a snow depth of ~ 0.5 m). The scene denoted 
by Derby Peak (39°59'N, 107° 10' W, ~ 3500ma.s.l.) is an- 
other high-altitude, rather flat region (in fact belonging to 
the Flat Tops Wilderness Area). The closest SNOTEL site, at 
Trapper Lake, 2950ma.s.l., reported a snow depth of 0.7 m. 


3 Methodology 

The method employed for the atmospheric correction pro- 
cedure is explained in detail in 012, to which the reader can 
refer if interested in the rigorous mathematical formalism. To 
provide a brief summary, the signal measured by the instru- 
ment is a nonlinear function of the surface reflectance due to 
the interaction of the radiation reflected by the surface with 
the atmosphere above. For this reason, a linear regression 
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cannot be applied to infer the model parameters. Neverthe- 
less, a cost function with the departure of the model from 
the measurements can still be computed and minimized, pro- 
vided the availability of the derivatives (Jacobians) of the ra- 
diative transfer model with respect to the descriptive param- 
eters. The Jacobians, generated without extra-computational 
cost from the linearization of the doubling-adding radiative 
transfer code in use at GISS (De Haan et al., 1987), are used 
in an iterative scheme based on an inversion and optimization 
procedure of the Gaussian-Newton kind (Rodgers, 2000) to 
search for a minimum of the cost function. The result is the 
decoupling of the total signal measured at instrument alti- 
tude into its surface and atmospheric contributions, i.e., an 
automated atmospheric correction. Within this scheme, the 
surface total reflectance is expressed as the weighted sum of 
three kernels as in the RossThick-LiSparse form (Wanner 
et al., 1995; Lucht et al., 2000): 


R(&) = /iso + /vol^volC©) + /geo^geoO), (5) 


where 0 is the scattering angle. The three terms on the right- 
hand side describe different scattering mechanisms: a Lam- 
bertian component for isotropic reflection, a volumetric term 
K v oi for scattering from small structures and a geometric 
term Kg eo for a macroscopic shadowing from large objects. 
Although this model was created for reflection from surfaces 
other than snow, its mathematical form adapts without prob- 
lems to snowpacks, with the isotropic kernel systematically 
assuming the largest weight because of the dominance of 
multiple scattering in the reflection of light from snow. 

As far as the polarized reflectance is concerned, we instead 
employ a scaled Fresnel kernel elsewhere used for vegetated 
surfaces (Nadal and Breon, 1999; Waquet et al., 2009): 


, F p (9j) 
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with the fractional polarized reflectance F p (0j ) given 
by (Born and Wolf, 1999): 
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where in our case n i and «2 are the indices of refraction of 
pure air and ice, respectively, and the angle of refraction 0 t is 
connected to the angle of incidence 0j by Snell’s law. 

For the retrieval of crystal habit, the polarized reflectance 
at 864 nm is used because of the very weak ice absorption 
at this wavelength and the relatively weak Rayleigh scat- 
tering contribution. Ice crystals are found to vary within a 
virtually countless set of shapes, but, limited by the lack of 
constraints on the correct optical properties, very few studies 
have departed from the common practice of assuming spher- 
ical shapes for the snow grains, an approximation which can 
negatively impact the quality of the albedo products (Tedesco 
and Kokhanovsky, 2007; Aoki et al., 2000). Among the avail- 
able optical models for non-spherical grains are for example 


those of Baum et al. (2011, 2014) and Yang et al. (2013). 
Incidentally, Raisanen et al. (2015) recently attempted to fit 
measured phase functions of blowing snow by using phase 
functions of several complex habits provided by Yang et al. 
(2013) in addition to phase functions of other habits. How- 
ever, the selection of habits in their study remained lim- 
ited, and to fit the data a rather arbitrary mixture of drox- 
tals, aggregates of plates and Koch fractals was needed, lead- 
ing them to conclude that this mixture “most probably does 
not represent properly the actual distribution of snow grain 
shapes in blowing snow (or snow on ground)”. In our case, 
the description of the snowpack is refined by treating it as 
a compact, optically semi-infinite medium composed of ice 
crystals modeled as hexagonal prisms with variable aspect 
ratios and microscale facet roughness. 

It should be noted that for the specific task of remote sens- 
ing of aerosols over snow, a model capable of accurately re- 
producing the contribution of the surface to the total signal 
should be considered successful regardless of its capability 
to mimic the true shape of the crystals. Our approach is based 
on the recognition that hexagonal plates and columns can be 
effectively used as radiative proxies for more complex crystal 
habits (Fu, 2007; Um and McFarquhar, 2007, 2009; Baran, 
2009; van Diedenhoven et al., 2012; Ottaviani et al., 2012a). 
As demonstrated by van Diedenhoven et al. (2012), match- 
ing measured polarized reflectances with a look-up table of 
values simulated by adjusting the aspect ratios and rough- 
ness parameters of simple hexagonal particles yields an esti- 
mate of the asymmetry parameter of more complex ice crys- 
tals. Our database, initially created to be of use for the re- 
trieval of ice cloud properties, contains 765 different combi- 
nations of aspect ratios and roughness parameters and was 
computationally assembled by running Monte Carlo simu- 
lations based on geometric optics (Macke et al., 1996), the 
performance of which was evaluated by van Diedenhoven 
et al. (2012). The retrieved aspect ratios are interpreted as the 
mean of the aspect ratios of the components of the complex 
crystals. 

The microscale surface roughness is statistically ac- 
counted for by tilting, for each interaction with a ray, the nor- 
mal to the crystal surface by an angle randomly selected in 
the S ■ [0°,90°] interval, where S is referred to as the rough- 
ness parameter. In our case, S is limited to 0.7 (i.e., a maxi- 
mum tilt angle for the facet of 63°) because for higher val- 
ues the probability of unphysical scattering events strongly 
increases, resulting in progressively larger loss of accuracy 
of the GO calculations of a given number of rays. Neshyba 
et al. (2013) have shown that various other definitions and 
parameterizations of crystal surface roughness lead to sim- 
ilar results. Yang et al. (2008) found that such approaches 
are efficient yet accurate treatments of microscale roughness. 
Furthermore, the effects of microscale surface roughness and 
macroscale crystal distortion were shown to be largely equiv- 
alent by Liu et al. (2014). 
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Figure 1. Atmospherically corrected polarized (second row) and total (third row) reflectance for three snowfields (columns) overflown in 
Colorado, USA. The first row shows real-time imagery from the high-resolution camera onboard the ER-2, captioned by the scene altitude 
above sea level together with the solar zenith angle and the Sun-sensor relative azimuth at the time of observation. The last row pertains to 
the search of the optimal ice grain model: the contour plots map the RMSE of the fit to the polarized reflectance at 864 nm when simulating 
the snowpack with each crystal habit (i.e., each combination of aspect ratio and roughness parameter) in the database. The red circles locate 
the values of the RMSE falling under a predetermined threshold and are taken to represent the retrieved (optimal) habits. 
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The retrieval of grain size is instead based on the total re- 
flectance at selected shortwave infrared (SWIR) bands where 
ice absorption is high. The total reflectance of an optically 
semi -infinite snow layer at 2264 and 1594 nm is effectively 
determined by the asymmetry parameter and the single scat- ^ 

tering albedo at the respective wavelength, which in turn is w 

mainly determined by the effective diameter of the ice crys- ^ 

tals. The grain diameter is not unequivocally defined in the < 

literature, which is understandable since the definition of < 

snow grain itself is often ambiguous. There is a consensus M 

for expressing grain effective diameter as proportional to the yj 

ratio of crystal volume to projected area, but the value of the ^ 

multiplier can vary according to the specific choice of the W 

author (Kokhanovsky et al., 2011; Kokhanovsky and Zege, O 

2004; Aoki et al., 2000). In line with several recent publica- 
tions (see, for example, Zege et al., 2011; Jin et al., 2008), 
we define it as 3/2 times the ratio of the crystal volume 
over its surface area, as outlined in Appendix A, which in- 
cludes the details of the relationship to geometric size. The 
effective diameter of snow for a given scene is estimated 
by matching the measured reflectance with an ice crystal 
model whose asymmetry parameter is consistent with that re- 
trieved for that scene using the method described above. To 
this end, various modified gamma size distributions were ap- 
plied to yield varying effective diameters, and the radiative 
transfer code was used to simulate the corresponding 1594 
and 2264 nm reflectances (van Diedenhoven et al., 2014). For 
conditions equivalent to those found in snowpacks, Bi et al. 

(2014) showed that the errors in retrieved effective radius at- 
tributable to the use of the conventional geometric optics ap- 
proach are below 5 %. Another important aspect to consider 
is that the 1594 and 2264 nm channels experience different 
penetration depths, with the former weighted more towards 
the top layer, while the latter probes deeper into the medium; 
as a consequence, the retrieved sizes contain information on 
the vertical structure of the snowpack (Wang et al., 2009; Li 
et al., 2001; Warren, 1982). 

As far as the atmospheric component is concerned, the 
high elevation of these alpine regions allowed us to assume 
negligible contributions from aerosols, so that in the radia- 
tive transfer code a standard Rayleigh atmosphere was used 
with a surface pressure extrapolated at the indicated altitude. 
Cirrus-free conditions were guaranteed from the screening 
procedure described earlier. The SWIR channels were cor- 
rected for the minimal interference caused by absorption 
from a standard amount of trace gases. 

4 Results and discussion 

The core of our findings is condensed in Fig. 1. For each of 
the three scenes (different columns in the figure) the analy- 
sis was performed over a single RSP pixel, corresponding to 
a spatial resolution of about 225 m. Indeed, averaging over 
a few adjacent aggregated scans did not show appreciable 



0.70 0.72 0.74 0.76 0.78 0.80 0.82 0.84 0.86 0.88 0.90 0.92 0.94 


Asymmetry Parameter, g(865) 



0.10 1.00 10.00 
ASPECT RATIO 


Figure 2. The “butterfly” pattern of the contour plot for the asym- 
metry parameter g as a function of aspect ratio and roughness pa- 
rameter resembles that obtained for the RMSE plots in Fig. 1. The 
red circles in Fig. 1 (here in green for Grand Mesa, magenta for 
Lake Granby and blue for Derby Peak) that identify the retrieved 
optimal crystal habits are found for each scene clustered around 
a single value of g. 


scan-to-scan variability, revealing more or less uniform snow 
conditions on a scale larger than the instrument instantaneous 
field of view. Real-time imagery from the high-resolution 
camera onboard the ER-2, overlaid to Google Earth, is in- 
cluded for context. The second and third rows of panels show 
the results of the iterative correction scheme applied to the 
relevant RSP channels. With the exception of 012, the visi- 
ble and SWIR behavior of the polarized reflectance of snow 
has not been previously published. 

The total reflectance exhibits the familiar high spectral 
albedo in the visible, albeit with some expected scene-to- 
scene variability. Lake Granby in particular shows lower to- 
tal reflectance than the other two scenes; in situ simultaneous 
surveys from the personnel at the US Forest Service in Ara- 
paho National Forest reported a snow depth on the order of 
5-8 cm. These values are close to the limit of thickness re- 
quired to visually mask the underlying surface, leaving the 
possibility that the lake ice underneath the snow cover is 
darkening the reflectance. Nevertheless, a layer of a few mil- 
limeters is sufficient to fulfill the semi-infinite approximation 
at the longer wavelengths used for the retrieval of grain size: 
the RSP SWIR channels at 1594 and 2264 nm are very dark 
because of the strong absorption of ice at these wavelengths. 

The polarized reflectance as a function of the scattering 
angle is very similar among the three scenes. Note that in the 
case of Derby Peak the flight track was oriented at about 45° 
with respect to the principal plane (i.e., the plane defined by 
the solar azimuth and the normal to the surface), leading to 
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a smaller range of scattering angles collected by the RSP. It 
has already been observed in 012 as the land model used for 
the polarized reflectance converges to zero toward backscat- 
ter (scattering angle > 160°), and in this region agreement 
with the data is not to be expected. It is useful to quantify 
the spectral spread of the channels by examining the differ- 
ences from the reference 2264 nm channel. The signal in this 
band is in fact a proxy for the surface contribution, since 
the Rayleigh (and fine-mode aerosol) contribution becomes 
more and more negligible as the observation windows shifts 
to the SWIR. If the spectral differences are small, that same 
signal can then be subtracted from the shorter RSP wave- 
lengths when aerosol retrievals are attempted. The residuals 
with respect to the 2264 nm polarized reflectance are found to 
be within 0.005, which as in 012 is larger than conventional 
land surfaces but still very small. The polarized reflectance is 
small too in absolute value. 

Shifting the attention to the search of a snow model, 
the plots at the bottom contour the root mean square error 
(RMSE) of the data fit to each combination of aspect ra- 
tio and roughness parameter in the database. Multiple com- 
binations led to a satisfactory fit to each of the three po- 
larized reflectance observations, so we considered as opti- 
mal all habits with an RMSE falling below a small thresh- 
old value (4.5 x 10 -7 ). These minima, marked with red cir- 
cles, therefore correspond to those crystal habits that when 
used to model the snowpack at the bottom of the simu- 
lated atmosphere yield the best fits to the data. The analysis 
shows that the polarized reflectances are consistent with par- 
ticles that have components that are like thin plates or like 
thin columns, a fact already observed in 012. Interestingly, 
the butterfly-like topology of the RMSE was recognized to 
closely resemble that obtained for the asymmetry parameter 
g plotted versus the same coordinates (Fig. 2). In each scene, 
the optimal crystal habits are associated with very similar 
values of the asymmetry parameter (0.84 for Lake Granby, 
0.876 for Derby Peak and 0.90 for Grand Mesa). This result 
is understood by considering that the asymmetry parameter is 
the main descriptor for the scattering properties. As demon- 
strated by van Diedenhoven et al. (2012) using a large variety 
of complex ice crystal shapes, all the combinations of aspect 
ratio and roughness parameter fitting a polarized reflectance 
measurement are ultimately characterized by similar values 
of g. Consequently, they will also yield similar particle size 
retrievals. 

Figure 3 shows the phase function and the degree of linear 
polarization, obtained for each scene by averaging the phase 
matrices associated with the retrieved minima. The lack of 
halo peaks in the 20-40° range, typical of the phase func- 
tions of pristine crystals, and the general smooth behavior 
(except from the strong forward peak) are consistent with the 
moderate-to-high roughness parameters determined in the 
analysis. 



Figure 3. Phase functions (upper panel) and degree of linear po- 
larization (DoLP, lower panel), resulting from the average over the 
optimal habits retrieved for each scene. The color coding is the same 
as for the circles in Fig. 2. 

All retrieved parameters are reported for comparison in Ta- 
ble 1, together with the results of the size retrievals. Using 
the simple approach outlined in the previous section, we ob- 
tained effective diameters equal to 152, 182 and 144 pm for 
the Grand Mesa, Lake Granby and Derby Peak scene, respec- 
tively, when using the RSP band at 2264 nm. Using instead 
the 1594 nm channel, diameters typically 30 pm smaller are 
obtained for each scene, which can be explained by the 
different penetration depths of the two channels, indicating 
crystal size increasing with depth as expected (Hori et al., 
2007; Li et al., 2001; Aoki et al., 2000; Warren, 1982). These 
values, at the lower end of common retrieval ranges, are nor- 
mally associated with fresh snow conditions as found in stud- 
ies based on other remote sensors in alpine regions (Negi and 
Kokhanovsky, 2011; Painter et al., 2009; Dozier and Painter, 
2004; Nolin and Dozier, 2000) and polar plateaus (Lyapustin 
et al., 2009; Hori et al., 2007). However, a direct compar- 
ison is challenged by differences in the definition of grain 
size, assumptions on grain sphericity and different penetra- 
tion depths achieved by the use of different instrument chan- 
nels. The larger grain sizes obtained for the Lake Granby 
scene were attributed to the shallow aged snow over the lake, 
since metamorphosed snow is recognized to be composed of 
larger grains (Aoki et al., 2000; Warren, 1982). 

5 Conclusions 

We have presented the extension of a novel analysis de- 
scribed by Ottaviani et al. (2012b) to a data set acquired with 
the high-accuracy RSP airborne instrument overflying high- 
altitude snowfields. The data were processed within an auto- 
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Table 1. List of ice crystals’ parameters retrieved for the three analyzed scenes. The value in parentheses indicates the RSP band used in 
the retrieval. The effective diameter can be obtained using either the 2264 or the 1594 nm channel (see text), together with a crystal habit 
characterized by the asymmetry parameter retrieved for the corresponding scene. 


Retrieved parameter 

Grand Mesa 

Lake Granby 

Derby Peak 

Asymmetry parameter (864 nm) 

0.90 

0.84 

0.876 

Aspect ratio (864 nm) 

<0.08 

<0.05 

< 0.2, >10 

Roughness parameter (864 nm) 

0.3-0. 5 

0. 6-0.7 

0.2-0. 6 

Effective diameter (2264 nm) 

152 pm 

182 pm 

144 pm 

Effective diameter (1594 nm) 

122 pm 

152 pm 

114 pm 


mated retrieval scheme capable of isolating the surface con- 
tribution to the total signal. The spectral dependence of the 
retrieved polarized surface reflectance is slightly larger than 
for soil or vegetated surfaces, but nonetheless small relative 
to typical aerosol contributions (cf. Fig. 6 in 012). This fact 
has important implications for the construction of advanced 
algorithms which aim at the retrieval of the microphysical 
properties of aerosol layers located over ice and snowfields. 

Furthermore, we applied a retrieval procedure which is 
demonstrated to yield snow microphysical parameters of pri- 
mary climatological interest. We retrieve snowpacks behav- 


ing as a collection of grains of extreme geometries (thin 
plates and/or long columns) and moderate to high microscale 
roughness. These results reinforce those of our previous 
study based on a limited data set acquired over very rugged 
terrain. For each scene, the best fits are obtained for a set 
of crystal habits characterized by very similar values of the 
asymmetry parameter and extreme aspect ratios. Size re- 
trievals were also tested, leading to grain effective diameters 
in the 140-180 pm range at the top of the snowpack. 
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Appendix A: Size descriptors for hexagonal prisms 


The results of Macke’s geometric optics code (Macke et al., 
1996) were tabulated with the choice of the radius of the 
sphere with equal surface area, r sp h, as a size descriptor. This 
parameter is linked to the geometric size of hexagonal prisms 
and to the grain effective diameter by the following relations : 


Ahex — 4 7T r“ ph 


(surface area of sphere and hexagonal prism), 
A p = Ahex/4 = 7t r“ ph 

(projected surface area of hexagonal prism), 


Hiex 


4 A n 


3 V3 + 12 a 

(radius or side length of hexagonal prism), 


V = 3 V3< 


hex 


(volume of hexagonal prism), 
6V3qi 


, 3 V 

4 '' = 2X = 


V3 +4cy 

(effective diameter). 


Hiex 


(Al) 

(A2) 

(A3) 

(A4) 

(A5) 


where a is the aspect ratio, defined as the ratio of the cross- 
section width of the prism to its length. The last expression 
is valid for a single particle of any solid shape. In case of 
a polydisperse collection of particles, the effective diameter 
is obtained through integration over a size distribution: 


3 f 0 °°V(D) N(D) dD 
2 / 0 °° Ap(D) N (D) d D 


(A6) 


with N denoting the number density of crystals with maxi- 
mum dimension D, which can be defined as for plates 
and 2 a r^ex for columns. 
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